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Abstract. This paper generahzes Moody's correlated binomial default distribution 
for homogeneous (exchangeable) credit portfolio, which is introduced by Witt, to the 
case of inhomogeneous portfolios. As inhomogeneous portfolios, we consider two cases. 
In the first case, we treat a portfolio whose assets have uniform default correlation 
and non-uniform default probabilities. We obtain the default probability distribution 
and study the effect of the inhomogeneity on it. The second case corresponds to 
a portfolio with inhomogeneous default correlation. Assets are categorized in several 
different sectors and the inter-sector and intra-sector correlations are not the same. We 
construct the joint default probabilities and obtain the default probability distribution. 
We show that as the number of assets in each sector decreases, inter-sector correlation 
becomes more important than intra-sector correlation. We study the maximum values 
of the inter-sector default correlation. Our generalization method can be applied to 
any correlated binomial default distribution model which has explicit relations to the 
conditional default probabilities or conditional default correlations, e.g. Credit Risk"*", 
implied default distributions. We also compare some popular CDO pricing models 
from the viewpoint of the range of the implied tranche correlation. 
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1. Introduction 

The modeling of portfolio credit risk and default correlation are hot topics and pose 
entirely new problems [H [21 [3], HI [5]. CDOs are financial innovations to securitize 
portfolios of defaultable assets. Many probabilistic models have been studied in order to 
price CDO tranches [6l[71[Hl[9l[ini[IIl[l2l[3l[T3l[Tl[T5]. Most of them are implemented 
with Monte Carlo simulations and as the number of names in a portfolio increases, 
the computational time increases. The Factor approach uses a small number of latent 
factors that induce the default dependency [11]. Conditionally on the latent variables 
values, default events are independent. It becomes easy to calculate the loss (default) 
distribution function. Along this line, some semi-explicit expressions of most relevant 
quantities were obtained [T6] . 

On the other hand, correlated binomial models were also proposed to describe the 
default dependency structures. The first one is a one-factor exchangeable version of 
CreditRisk [TTl [181 [I9]. The aggregate loss distribution function is given by the beta- 
binomial distribution (BBD). The second one is Moody's correlated binomial default 
distribution model, which was introduced by Witt [20]. ( hereafter the MCB model) 
The authors also consider the applicability of the long-range Ising model [211 [22]. These 
models use Bernoulli random variables. Differences stem from different definitions of the 
conditional correlations [23]. In the MCB model [20], the conditional default correlation 
between assets is set to be constant irrespective of the number of defaults. Those of 
BBD decay with an increase in default. We are able to adapt a suitable form for the 
conditional correlations. Recently it has become possible to calibrate pn from implied 
default distributions [2l]. By using the "implied correlated" binomial model, whose 
conditional correlations are those of the implied distribution, it may become easy to 
estimate hedge ratios and so forth. 

The advantage of these correlated binomial models come from the fact that they 
are easier to evaluate than other more refined models. If a probabilistic model is 
implemented by a Monte Carlo simulation, the evaluation of the price of these derivatives 
consumes much computer time and the inverse process to obtain the model parameters 
becomes tedious work. With the above correlated binomial models, one can estimate 
the model parameters from the CDO premiums more easily. However, these models are 
formulated only for homogeneous portfolios, where the assets are exchangeable and they 
have the same default probability p and default correlation p. Generalization to more 
realistic inhomogeneous portfolios where assets have different default probabilities and 
different default correlations should be done. 

In this paper, we show how to generalize Moody's correlated binomial default 
distribution (MCB) model to two types of inhomogeneous portfolios. Our generalization 
method can be applied to other correlated binomial models, including implied correlated 
binomial distributions, by changing the condition on p„. We obtain the default 
probability function PN{n) and examine the dependence of the expected loss rates of the 
tranches on the inhomogeneities. With the proposed model, we also estimate the implied 
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values of the default correlation. Comparison of the range of the tranche correlations 
with those of the Gaussian copula model and BBD model are also performed. 

About the organization in this paper, we start with a short review of Moody's 
correlated binomial default distribution (MCB) model in Section [2l The dependence 
of pn on the number of defaults n is compared with the BBD and Gaussian copula. 
Section [3] is the main part of the paper. We show how to couple multiple MCB models 
as a portfolios credit risk model for inhomogeneous portfolio. In the first subsection, we 
couple two Bernoulli random variables X, Y and recall on the limit of the correlation 
Pxy between them. In the next subsection, we couple an MCB model of assets 
Xi,X2, - ■ ■ ,Xn with a random variable Y and we study the maximum value of the 
correlation between Xi and Y. Then we couple two MCB models with and M 
assets. Choosing the model parameters properly, we construct an MCB model for an 
inhomogeneous default probability case and obtain the default distribution function. 
The last subsection is devoted to an inhomogeneous default correlation case. Assets 
are categorized in different sectors and inter-sector and intra-sector default correlations 
are not the same. We consider a portfolio with K sectors and k—th sector contains A^fc 
assets. Within each sector, the portfolio is homogeneous and it has parameters as pk and 
Pk- The inter-sector default correlations are not the same and they depend on the choice 
of sector pairs. We construct the joint default probabilities and the default probability 
function P/v(?2) for the portfolio explicitly. In Section HI using the above results, we 
estimate the implied default correlation for each tranche from a CDO's market quotes ( 
iTraxx-CJ Series 2). We compare the range of the correlations of MCBs, BBD and the 
Gaussian copula model. We conclude with some remarks and future problems. 

2. Moody's correlated binomial default distribution 




Figure 1. Homogeneous portfolio with N assets. The assets are exchangeable and 
the default probability is p and the default correlation is p. The state of i-th asset is 
described by a Bernoulli random variable Xi. 

We review the definitions and some properties of Moody's correlated binomial 
default distribution (MCB) model. We consider a homogeneous portfolio, which 
is composed of exchangeable A^ assets. Here the term "homogeneous" means that 
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the constituent assets are exchangeable and their default probabilities and default 
correlations are uniform. We denote them as p and p. Bernoulli random variables 
Xi show the states of the i-th assets. Xj = 1 means that the asset is defaulted and the 
non-default state is represented as Xi = 0. The joint default probabilities are denoted 

as 

P{xi, X2,--- , xn) = Prob(Xi = xi,X2 = X2, - ■ ■ ,Xn = xn)- (1) 

In order to determine P{xi,X2, ■ ■ ■ ,xn), we need 2^ — 1 conditions for them. Here 
2^ corresponds to the number of possible configurations and —1 comes from the 
overall normalization condition for the joint probabilities. From the assumption of the 
homogeneity for the portfolio, the number of degrees of freedom of the joint probabilities 
are reduced. The probability for n defaults and N — n non-defaults is the same for any 
configuration (xi, X2, ■ ■ ■ , xat) with Y.f=iXi = n. The number of defaults n is ranged 
from to and considering the overall normalization condition, remaining degrees of 
freedom are A^. 

In the MCB model, the conditional default probabilities are introduced. We denote 
Pn as the default probability for any assets under the condition that any other n assets 
of the portfolio are defaulted. To exemplify the situation concretely, we take the n assets 
as the first n oi N assets, we denote them with n' as n = 1, 2, ■ ■ ■ , n. The condition 
that they are defaulted is written concisely as O^'^i^n' — 1- The conditional default 
probability for n + 1-th assets under the condition of n defaults can then be written as 

n 

Pn=<Xr,+l\ n^„'=l>- (2) 

n'=l 

Here, < A\C > means the expected value of random variable A under condition C 
is satisfied. Xn takes 1 for n-th asset default, < Xn\C > corresponds to its default 
probability under condition C. Of course, any asset from k = n + 1, ■ ■ ■ ,N can be chosen 
in the evaluation of the expected value for p„ under the condition that D^'^i -^„' = 1- 
Po is nothing but the default probability p. 

N independent conditional default probabilities Pn{i = 0, ■ ■ ■ , N — 1) are determined 
by the following condition on the default correlations. 

Cor(X,+i,X„+2| n X^> = l) = p. (3) 



n =1 



Here, Cor(X, F|C) is defined as 



Cor(X, Y\C) = <XY\C>-<X\C><YIC> 

'< X\C > (1- < X\C >) < Y\C > (1- < Y\C >)) 



The conditions on the default correlations give us the following recursion relations for 

Pn as 

Pn+l = Pn + - Pn)p- (5) 

These recursion relations can be solved to give Pn as 

p„ = l-(l-p)(l-p)". (6) 
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Pn increases with n and p„ — ^ 1 as n — oo for p > 0. 

From these conditional default probabilities Pn{n = 0, ■ ■ ■ , N — 1), the joint default 
probabilities for the configuration x = {xi,X2, ■ ■ ■ ,xn) are given as 

N 

Fix) = Pix^,X2, ■ ■ ■ =< n ^^(1 - ^n)'-"" > . (7) 

n=l 

The normalization condition for P{x) is guaranteed by the following decomposition of 
unity. 

N N 1 N 

1 =< 1 >=< n{x.+(i-xo} >= n [ E ] < n ^^-ii-^j^-^-- > (s; 

71=1 j = l Xj=0 71 = 1 

The probability for n defaults is 

n N 

ptvw =7vc„ X p(i,i,---,i,o,---,o) =^c„ < n (l-^^)> 

i=l 4 = 71+1 

N-n 71+fc-l 

= NCnY.N-nCk{-l)\ H Pn') ■ (9) 
fc=0 n'=0 

We modify the above MCB model as follows. In the MCB model, the default 
correlation is set to be constant irrespective of the number of default (see([3])). We 
change the condition as 

71 

Cor(X„+i,X„+2| n ^n' = 1) = pexp(-nA). (10) 

n'=l 

Here, we introduce a parameter A > and the default correlation under n defaults decay 
as exp(— nA). If we set A = 0, the modified model reduces to the original MCB model. 

There are two motivations for the modification. The first one is that it is 
mathematically necessary to couple multiple MCB models. We discuss the mechanism 
in greater detail in the next section. Here, we only comment on the limit value of pn as 
n ^ oo. The modification changes the recursive relation for p„ to 

Pn+l = Pn + {I - Pn)pexp{-nX). (11) 

Pn is calculated as 

n-l 

P„ = l-(l-p) n(l-Pn')- (12) 

n'=0 

Here p„ is defined as p„ = pexp(— riA). 

Pn increases with n, however the increase is reduced by the decay of the correlation 
with n. The limit value of Pn with n ^ cxd is roughly estimated as 

Poo = lim p„ = 1 - (1 -j9)(l -p)exp(-p- -). (13) 

n — ^oo — 6~ 

For A = 0, Poo = 1 and Poo = P + (1 — p)p = Pi for A = oo. In Figure [2|, we show the 
enumerated data for poo and the results from eq. ffTSl) . As A increases, Poo decreases and 
the A dependence is well described by eq.( !T3|l . 
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Coditional Default Probability Pp at n=lnfinity 
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The second motivation is that popular CDO pricing models have decaying 
correlation with n. The BBD model's p„ is given as 



Pn = (14) 

1 + np 

In the Gaussian copula model, we do not have the explicit form for From its 
aggregate loss distribution function, it is possible to estimate them. In Figure [3l we 
show Pn for the MCB, BBD and Gaussian copula models. We set p = p = 0.1 and 
= 30. The Gaussian copula's p„ does not show monotonic dependence on n. After a 
small peak,, it decays to zero. In order to mimic the Gaussian copula model within the 
framework of a correlated binomial model, such a dependence should be incorporated 
in the assumption on p„. 

Figure H] depicts pn in the same setting. With the same p and p, all pn curves 
pass through pq = p at n = and pi = p + {1 — p)p at n = 1. After n > 1, the 
behaviors of p„ depend on the models' definitions on p„. p„ saturate to about 0.4 for 
MCB(A = 0.3), which means that a large scale avalanche does not occur and the loss 
distribution function has a short tail. In MCB with A = 0.0 and Gaussian copula 
models, their p„ saturate to 1. The behaviors are reflected in the fat and long tails in 
their loss distribution. 

Figure [5] shows the semi-log plot of P'iQ{n) for the MCB, BBD and Gaussian copula 
models. We also plot the binomial distribution Bin(30,0.1). The default correlation 
shifts the peak of the binomial distribution to n = and Pzo{n) comes to have a long 
tail. MCB, BBD and Gaussian copula have almost the same bulk shape. In particular, 
in MCB, even if we change A, PzQ{n) has almost the same shape for n < 15. The 
bulk shape of PN{n) is mainly determined by p„ with small n. p„,s with large n comes 
from very rare events n"'=i ^n' ~ ^ ^^"^ contains information about the tails of the 
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Figure 4. p„ for MCB (solid, dotted), BBD (+) and Gaussian copula models (x). 
We set p = p = 0.1 and N = 50. 



distributions. They do not affect the bullc part significantly . 

There are differences in their tails. One sorts the models in the order of thinnest 
tail to fattest tail, we have 

MCB(A = 0.3) < BBD < Gaussian copula < MCB(A = 0.0). 

MCB(A = 0) has almost the same shape as the Gaussian copula. However, it has a 
bigger tail than the gaussian coupla at n = 30. The tail of MCB(A = 0.3) is short 
compared with other models. We can understand this behavior from the behavior of p„. 
We also note another role of the damping parameter A. In the calculation P/v(n), 
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Figure 5. Loss Distribution P3o{n) for p = 0.1 and p = 0.1. A = 0.0,0.3 and 
p = p = O.l. We also plot the profile for the binomial distribution. 



there are many cancellations Y.k=o' N-n.Cki—^)'' in (ED from the decomposition of 
niln+i(l ~ -^i) ■ This causes numerical errors in the evaluation and it is difficult to 
get PAr(ri) for > 40, even if we use long double precision variables in the numerical 
implementation. When we set A = 0.3, the numerical error diminishes greatly and we 
can obtain P/v(n) even for = 100. This point is important when one uses the MCB 
model for analysis of the actual CDOs that have at least 50 assets. In addition, with 
A > 0, we can take p to be negatively large enough. In S&P's data, a negative default 
correlation of 0.1% or so has been reported [25]. We think that this point is also an 
advantage of the modified model. 

Hereafter, we mainly focus on the generalization of the MCB model. However, the 
same method and reasoning should be applicable to other correlated binomial models 
with any assumption on p„. If we set p„ as in (fT4|) . we have Beta-Binomial default 
distribution models for inhomogeneous portfolios. 

3. Generalization to Inhomogeneous Portfolios 

In this section we couple multiple MCB models and construct the joint default 
probabilities and P/v(n) for inhomogeneous portfolios. In addition, we show that the 
inter-sector default correlation can be set to be large enough by choosing A and other 
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parameters. We think that it is possible to use the model as a model for portfolio credit 
risk. 



3.1. Coupling of X and Y: 1+1 MCB model 

Before proceeding to the coupling of multiple MCB models, we recall some results for 
the coupling of two random variables X and Y . The default probability is Px iPy) foi' 
X (respectively for Y) and the default correlation between them is p^y, 

Corr(X, Y) = p^y. (15) 

As in the MCB model, we introduce the conditional default probabilities as 

Po = Px =< X > and pi =< X\Y = 1 > . 

Qo = Py =< y > and qi =< Y\X = 1 > . (16) 
From the default correlation pr^y, pi and qi are calculated as 



Pi=Po + (1 -Po). 



Po(l - go) 
\\ {l-po)qo 



Pxy 



gi = go + (1 - qo)> 



9o(l -Po) 
\{l-%)Po^''- 



(17) 



In the symmetric (homogeneous) case p^ = Py, the equality pi = qi holds and they are 
given as 

Pi = qi=Px + i^-Px)pxy 

The correlation p^y can be set to be 1 and in the limit pi = qi — > 1. The maximum 
value of Pscy is 1 in the symmetric case. Conversely in the asymmetric case {p^ ^ py), 
Pxy cannot set to be 1. The maximum value of pxy is determined by the condition that 
Pi < 1 and qi < I. From these conditions, we derive the following conditions p6] on 

Pxy as 



P: 



xy 



< 



go(l -Po) 



and Pxy < 



^ (1 - qo)po 
We introduce an asymmetric parameter r as 

r = Py 

Px 

and a function f{r,p) as 



Po(l - go) 
\ (1 -Po)go' 



f{r,p) 



(1 — p)r 



:i8) 



(19) 



(20) 



/ 1 — rm 

The maximum value of pxy is then given as 

Max(p,,) = Min(/(r,po),/(r,Po)"'). (21) 

Here Max{pxy) represents the maximum value of pxy and Min(A, B) means that the 
smaller value of A and B is taken. Figure E] shows Ma.x{pxy) as a function of the 
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Maximum value of Pj,^ vs r 
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Figure 6. Ma.x{pxy) for ~ 0.1(solid line), 0.5(dottcd line) and 0.2 < r < 5.0. 



asymmetric (inhomogeneity) parameter r. We show two curves, the sohd one for p^. = 0.1 
and the dotted one for px = 0.5. As the inhomogeneity r increases, that is r departs 
from r = 1, Maxlp^y) decreases. For fixed r, as px becomes large, Ma.x{pxy) becomes 
smalL The reason is that the condition pi < 1 becomes more difficult to satisfy as 
Px = Po increase, pi is a monotonous increasing function of po. In the previous section, 
the conditional default probability Pn becomes smaller as we set A larger. When we 
set a large A, we show that it is possible to couple multiple MCB models with strong 
default correlation. 



3.2. N + 1 MCB model 

For the second step, we couple an assets MCB model with one two-valued random 
variable Y. We introduce N random variables X„(n = 1,---,A^) and the default 
probability and the default correlation for them is < X„ >= Px and p. The default 
probability for F is < Y >= py and the default correlation between Y and Xn is 
written as pxy We assume homogeneity for the assets MCB model and the default 
correlation between X„ and Y is independent of the asset index n. As in the previous 
cases, we introduce conditional default probabilities as 

n n 

Pnfl =< Xn+i\ n ^n' = 1 > and pn,i =< A:„+i| Y[ X^' X Y = 1 > 

71 =1 71 =1 

71 

qn =< F| n ^n' = 1 > ■ (22) 

n'=l 
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The joint default probabilities P{x,y) = P{X = x,y) are calculated by decomposing 
the following expression with these conditional default probabilities 

N 

Fix, y) =< n ^^"(1 - ^n)'""" X - Yf-y > (23) 

71=1 

The joint default probabilities with the condition Y = 1 are 

N 

P(x\Y = 1) =< n X^-{1 - X^y--- > /py. (24) 

71=1 

The joint probabilities with the condition F = are obtained by the following relations. 

N 

P{x\Y = 0)(l-p,)+P(f|F = l)py = P{x) =< n X^"(l-Xji-^" > (25) 

n=l 

The explicit form for P{x\y = 0) is 

PfsiY = 0) = ni)-p(^ = DP,^ PS) 

1-Py 

The probability for n defaults is given as in eq.([n]) 

N-n n+k-1 

PN{n) = nCu E N-nCu{-l)\ n ^'-'.o)- (27) 

fc=0 n'=0 

The probability for n defaults with F = 1 is 

N-n n+k-1 
PN{n\Y = 1) = E N-nCk{-l)\ H Pn',l). (28) 

fc=0 n'=0 

Using the same argument as for the joint probabilities with F = 0, P]s[{n\Y = 0) is 
written as 

= 0) = (29) 

l-Py 

About the conditional default probabilities Pnfly we impose the same conditions as 
with the homogeneous assets MCB model. 



Cor(X„+i, X„+2| n ^n' = 1) = Pe-"" (30) 

n'=l 

The same recursive relation (ITT]) for p„ o is obtained and the Pnfl is given by 

n-l 

Pn,0 = I - {I - Px) Y[{1 - Pn')- (31) 
n'=0 

Here p„ is defined as before. Fot the conditions on p„ i and qn, there are two possible 
ways to realize a strong correlation. The only way to realize a strong correlation p^y 
between Xn+i and Y is 

n 

Corr(X„+i, Y\ H ^n') = Pxye~"\ (32) 

n'=l 

In this case, the relations for pn,i and g„ are 



Pn,iqn = qn+lPnfi = Pn,Oqn + PxysJPnfii^ -Pn,o)?n(l - qn)e (33) 



Moody's Correlated Binomial Default Distribution 
The recursive relations are 



\ (1 -Pn,o)qn 



qn+1 = qn + Pxye "^(1 - g„ 



(1 -Pn,0)qn 
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(34) 



(35) 



^ Pn,o(l - qn)' 

If we set Py = px and p^y = p, these relations reduce to pn,i = Pn+i,o and qn = Pn,o and 
this coupled model is nothing but the (A^ + 1) assets MCB model. 



Maximum value of p„ vs r 




Figure 7. Max{pxy)/p vs ; 
and the dotted line A = 0.3 



_ Py 



N ^30, 0.1, px = 0.1. The solid line A = 0.0 



We write the ratio as r„ and the conditions that pn,i < 1 and qn+i < 1 are 
summarized as 

pxy < e"^ X Min(/(r„,,p„,o), /(?^n,,Pn,o)~^)- (36) 
Mm{f{rn,Pn,o), f{rn,Pn,o)~^) is nothing but the condition for Max^p^y) of the two 
random variables X,Y with default probabilities {px,Py) = {Pn,o,i^nPn,o) (see eq. fl2T|) ). 
As explained above, Min(/(r„,p„), /(r„,p„)~^) takes maximum value 1 at r^i = 1 for 
any value of p„. It also decrease with the increase of p„ for fixed r„. 

The necessary condition for the model to be self-consistent is that Pn.i < 1 and 
qn < 1 for all n. We discuss A = and A > cases separately. 

• A = 0.0: Pnfl increases with n and Pn,o ^ 1 as n ^ oo. The range of Pnfl is 
[Px,Pn-i] — [Px; 1] and it is difficult to choose Py such that r„ ^ 1 for all n. If 
r„ departs much from 1 for some n, Ma.x{pxy) decreases. The choice px = Py and 
Pxy = p is possible, we anticipate that Max(p2;j^) decreases from p as Py departs 
from px . 
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• A > 0: The limit value of Pn becomes small (see eg. ( 1131) ) and the range of j9„ is 
narrow as compared with the A = 0.0 case. It may be possible to choose Py such 
that the asymmetric parameter r„ is small for small n. In addition, for large n, it is 
easy to satisfy the condition (136!) because of the prefactor exp(nA). We think that 
Max^pxy) is large in this case. 

We have checked numerically the values of the joint probabilities for all 
configurations {x,y) with iV = 30 and px = p = 0.1. In figure U\ we show the data 
of Ma.x{pxy) / p for the case A = 0.0 (solid line) and A = 0.3 (dotted line). In the A = 0.0 
case, Max{pxy) — p near r = 1 and as r departs from 1, Max{pxy) decreases from p. The 
data for A = 0.3 case show that it is possible to set a large Ma.x{pxy) if we use a large 
r. We can set p^y as strong as several times of p. 

We also point out that the above + 1 MCB model can be used to describe a 
credit portfolio where one obligor has great exposure. Such an obligor is described by 
Y and other obligors are by Xn- If Py is quite different from px, we can couple Y and 
Xn with strong p^y — Px by setting a sufficiently large A. 



3.3. N + M MCB Model : Coupled MCB model 




Next we consider a portfolio with two sectors. The first sector has assets 
and the second has M assets. To construct the joint default probabilities for the 
portfolio, we try to couple two MCB models. The former model's TV assets are described 
by Xn{n = 1,---,N) and the states of the latter model's assets are described by 
Ym{m = 1, ■ ■ ■ , M). The default probability and the default correlation in each sector 
are {px,Px) and {py,py). The default correlation between the assets in different sectors 
is denoted Figure [HD . 

Introducing the conditional default probabilities pn,m and gn,mas 

n m 

Pn,m = < Xn+l\ Y[ ^n' ]]_ Y^' = 1 >, Pofl = Px (37) 
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n m 

qn,m = < Ym+l\ ^n' ]]_ Ym' = I >, QOfi = Py, (38) 

n'=l m'=l 

we impose the following conditions on p„ o and go,m 

n 

Cor(X„+iX„+2| n = 1) = Pxexp{-n\x) (39) 

n'=l 
m 

Cor(ym+iFm+2| n Ym' = I) = Py(ixp{-m\y). (40) 

m'=l 

The recursive relations for pnfl and go,m are 

Pn+1,0 = Pn,0 + Px exp{-n\x){l - Pn,o) (41) 

go,m.+i = go,n + Pyexp{-mXy){l - go,m)- (42) 
Their solutions are, by denoting p^. ^ = p^. exp(— nA^.) and Py^m = PyGxp{—mXy), 

n-l 

Pn,0 = 1 - (1 - Px) n (1 ~ P2;,n') (43) 
n'=0 

go,m = l-il-Py) n (1 - P?/,-')- (44) 

m'=0 

For the inter-sector correlation, we impose the next conditions on Pn,m and qn,m, which 
is a natural generalization of + 1 case (see eq. (l32l) ). 

n m 

Cor(X„+iF„+i| n Xn' n = 1) = Px,e-("^^+'"^-). (45) 

n'=l m'=l 

We obtain the following recursive relations. 



Pn,m+lQn,m Qn+l,mPn. 



m 



Pn.mQn.m ~l~ Pxy^ ^ ^Jpn,m ( 1 Pn,m,)'?n,m(l Qn^m) • 

(46) 

Using these relations, we are able to calculate p„^m and g„^m iteratively starting from 
p„,o and go,m- 

The joint default probability for the portfolio configuration (x, y) is calculated by 
decomposing the following expression with Pn,m and 

P(f , y) = P(xi, X2, ■ ■ ■ , Xat, yi, 1/2, ■ ■ ■ , Z/m) 

= < n ^^"(1 - ^n)'-"" n ^m'"(l - ^™)'""" > • (47) 
n=l m=l 

In particular, the probability for n, m defaults in each sector, which is denoted as 
PN,M{n,m), is 

PN,M{n,m) = X mC„ X P(l,l,---,l,0,---,0, l,l,---,l,0,---,0) 

n N m M 

= NCmy<MCmy<< U^n' U (1 - X,) l[ Yr^' H ' ^l) > 
n'=l fc=n+l m'=l /=m+l 

N—n M-m 

= X uCm ^ X! X! (~^) ^ N-nCk X M-rnCl 

k=0 1=0 

n+k m+l 

X < n n > • (48) 
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PN+Niin) is easily calculated from Pp^^Min.m) as 

n 

PN+M{n) = Pn,nM'^ n-n'). (49) 

n'=0 

In decomposing < nn=i Ytm=i > , one can do it in any order. The independence 
of the order of the decomposition of < 11^=1 Ylm=i > is guaranteed by (145|) and 
( 146|) . We decompose it as 

k I fc-l l-l 

< Yl Xn Y[ Ym >= Y[ Pnfi X H ^k,m- (50) 

71=1 7Tt = l 71=0 7Tl=0 

For the maximum value of p^y, it is necessary to check all values of the joint 
probabilities. However, N + M model is reduced to + 1 or 1 + M model by choosing 
M = 1 or = 1 respectively. From the discussions and the results in the previous 
subsection for the + 1 model, we can anticipate as follows. 

(i) Xx = Xy = 0: If Px = Py and px = Py, two MCB models are the same and we can 
set pxy = Px = Py Two models merge completely and we have a (A^ + M) assets 
MCB model. As the asymmetry between the two models becomes large {px 7^ Py 
or Px 7^ Py) , Ma.x{pxy) decreases from Px, py 

(ii) Ax, Ay > 0: As Aa;, Aj; increase, Ma^yi^pxy) becomes large. The asymmetry in Px.Py 
and Px, Py diminishes yiax.{pxy)- 

Using the above coupled N + M MCB model, we study the effect of the dispersion 
of the default probability on Pis[{n) and on the evaluation of tranches. More complete 
analysis about the difference between the usage of individual spreads and of portfolio 
average spreads in CDO pricing has been performed in [5]. There, the usage of the 
average spread results in the lower estimation of the equity tranche. We consider a 
portfolio with N + N assets. The assets in each sector have default probabilities p± AP^ 
and an intra-sector default correlation p. We set the inter-sector default correlation pxy 
also as p. The inhomogeneity in the default probability is controlled by AP^. If we set 
AP^; = 0, the two sector are completely merged to one sector and we have a homogeneous 
2N MCB model with p and p. In figure [H we shows the default probability difference 
AP{n) between the inhomogeneous case P/v+Ar(n) with AP^ 7^ and homogeneous case 
P2Ar(n). AP(n) is defined as 

AP(n) = PN+Nin) - P2N{n). (51) 

We set A^ = 20, p = p = 0.03 and A = 0.3. The solid curve represents the data for 
APd = 0.01 and the dotted curve stands for the case AP^ = 0.02. We see that AP(n) is 
large only for small n. We also plot the results for the N + M BBD model. To construct 
the loss distribution function, it is necessary to change ( l39l) and ( HOl) to 

n 

Cor(X,+iX,+2| n ^n' = 1) = p/(l + np) (52) 

n'=l 
m 

C0l{Y^+^Yrn+2\ n = 1) = + ^P)- (53) 

m'=l 



Moody's Correlated Binomial Default Distribution 



16 



0.02 
0.015 
0.01 

0.005 



Q. 

-0.005-|t 



-0.01 



-0.015 -I 



-0.02 



-0.025 



X 



Effect of Dispersion AP^j on MCB and BBD 



IVICB APd=0.01 
IVICB APd=0.02 
BBD APjj=0.01 
BBD APcj=0.02 




— )K — )K — ) 



10 



12 



14 



Figure 9. AP(ri) vs n. Wc set N ^2Q,p^ 0.03, p ^ 0.03 and A^^ = A,, = 0.3. MCB 
in solid and dotted lines, and BBD in + and x symbols. 



Eq. (HSl) is also changed as 

n m 

Cor(X„+iF,,,+i| n Xn' n = 1) = + + ^)P)- 



(54) 



n'=l m'=l 

The recursive relations are also changed, but the remaining procedures are the same as 
for the N + M MCB model. As we have shown in the previous section, the bulk shapes 
of the loss distribution of MCB and BBD models are almost the same, and the effects 
of APd on them are also similar. 

To see the effect of AP^ on the evaluations of tranches, we need to study the change 
in the cumulative distribution functions D(i). D(i) is defined as 



N 



(55) 



D{i) represents the expected loss rate of the i—th tranche, which is damaged if more 
than i assets default [22] . One of the important properties of D{i) is 

1 N 



N ' 



p. 



(56) 



i=l 



This identity means that the tranches distribute the portfolio credit risk between them. 
Another important property of D[i) is that the expected loss rate D[i,j] of the layer 
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AD(i) and AD(i)/D(i) vs i 



Q 
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Figure 10. Plots of AD{i)/D{i) vs APd and AD{i) vs APd. N = 20,p = p = 0.03 
and APd = 0.01. We plot for 1 < n < 15. 



protection with attachment point i and detachment point j can be built up from D{i) 
as 



(57) 



3 

Figure [10] shows the plots of /AD{i) and of the ratio AD{i)/D{i) vs i. We see 
that the lower tranches' (z < 2) expected losses increase by APd, which is reasonable. 
The default probability of half the assets of the portfolio increases p ^ p + AP^ and 
the expected losses of the subordinated tranches increase. For the senior tranches, the 
absolute value of AD{i) decreases with i, however this does not mean that APd does 
not have a small effect on them. The absolute value of D{i) also decreases with i. From 
Figure fTOl we see that the ratio AD{i)/D{i) does not necessarily decrease with i. 



3.4- Multi Sector Case 

To couple two or more MCB models, we consider a portfolio with assets, which are 
categorized in K different sectors. Figure fTTl sketches the structure of the portfolio. The 
/c-th sector contains assets and the relation X^fcLi = N holds. The states of the 
assets in the /c-th sector is described by X!^^{nk = 1, ■ ■ ■ , Nk) and the default rate and 
default correlation are denoted as pk and pk- For the inter-sector default correlation, we 
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denote this as pij for the default correlation between the i— th and j— th sector. The 
intra-sector default correlation and inter-sector default correlation are different and the 
former is larger than the latter in general |27j . 




Figure 11. Structure of the portfolio. There are K sectors and fc-th sector contains 
Nk assets. The default probability and default correlation for the assets in the fc-th 
sector is pk, Pk- The inter-sector default correlation between the i— th and j-th sector 
is Pij. 

We have not yet succeeded in the coupling of three or more MCB models by 
generalizing the result for the coupled N + M MCB model. The reason is that the 
self-consistency relations are very rigid restrictions on the MCB model. These relations, 
for which the two-sectors version is given by (H6ll . assure the independence of the 
order of the decomposition in the estimation of the expected values of the product 
of random variables. It is difficult to impose any simple relations on the conditional 
default probabilities that satisfy the self-consistency relations. 

In order to construct the joint default probabilities for the assets states 
[x^, x^, ■ ■ ■ , x^), we do not glue together K MCB models directly. Instead, as depicted 
in figure [121 we glue multiple MCB models through one random variable Y. More 
concretely, we prepare K sets of A'^fc + 1 MCB models. A'^^ + 1 MCB model is the 
A^fc MCB model coupled with Y. The probability of P{Y = 1) is written as Py. We 
Introduce the following conditional default probabihties 

pLo =< n = 1 > p^^,, =< x„v,i n K', xY=i> 
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4 =< y\ n K', = 1 > • (58) 

We also impose the following conditions on Pn,.fiyPnk,i Ink 

Cor(X„\+„ n = 1) = Pfc exp(-n,A) (59) 



Corr(X„\+„y| n X;,) =p,,exp(-n,A). (60) 



The joint default probabilities P{x^, y) and the conditional joint default probabilities 
P^{x^\y) are constructed as before. 

P(x^ y) =< n «)^- (1 - X F^(l - Yy~y > (61) 

p(x^y) = p'=(x'=ll/)p(y). (62) 

Packing these conditional default probabilities P''{x'^\y) into a bundle, we construct the 
joint default probabilities for the total portfolio as 

p(x^f^■■■,f^)= 5: p(y)x n^(44---,4jy)- (63) 

y=Q,l k=l 

We also obtain the default probability function P/v(ni, 77-2, ■■■ , uk) for rifc default in the 
k-th sector as 

K 

PN{ni,n2, - ■ ■ ^tik) = P(y) ^ I[PNkink\y) (64) 

y=0,l k=l 

From the expression, it it easy to calculate the probability for n defaults and we write 
it as Pj^in). 

For the default correlation between the different sectors, we can show the next 
relations. 

Pij = Piy X pjy. (65) 

More generally, the conditional inter-sector default correlations obey the following 
relations. 

Cor(X;+„X4.+,| n K'^ n K' = 1) = P.,exp(-(n, + n,)A) (66) 

n'=l ri'=l ^ 

These relations mean that our construction procedure is natural from the viewpoint 
of the original MCB model. In particular, in the K = 2 case, these relations are 
completely equivalent with those of the N + M coupled MCB model. See (1451) and fl66|) . 
The conditional default probabilities obey the same conditions. 

Furthermore, from the results on Max{prcy) of the + 1 MCB model, we see that 
the model can induce a realistic magnitude of the inter-sector default correlation. By 
choosing A and Py properly, it is possible to set pky as large as several times of pk- 
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Figure 12. Gluing multiple MCB models with Y. The correlation between X^^^ and 
Y is Pky ■ The default correlation between X^, and X^^^ is given as pij — piy x pjy . 



Figure [13] plots Max(pa;y)=Max(pfcj/) as functions of p = pk- We set parameters as 
depicted in the figure. The sohd hne depicts the data for f = ^ = 3.0. The other two 
curves correspond to r = 1.0 (+) and r = 15.0 (x). We see that by setting r = 3.0. 
it is possible to set pky as large as several times of pk- If the intra-sector correlation 
Pi = pj is 10%, we can set piy = pjy = 20%. The inter-sector correlations is then 
Pij = Piy X Pjy = 4%. In general, pinter is smaller that Pintra, we think that the present 
model can incorporate a strong enough inter-sector default correlation. 

To prove the relations (!65l) and (!66|) . in order to calculate the correlation, we need 
to estimate the next expression. 

< Xi,^^,Xi^^,\ l[Xl,llX^, = l>. (67) 

i'=i j'=i 

If we fix the random variable Y, X^. and are independent. They are coupled by Y 
and the above equation is estimated by the average over Y = and y = 1 as 

< ^n,+i^ij+i\ n -^i' n -^j' = 1 > 

i'=i j'=i 

= < ^n.+ll n = l,Y = l>< X^.+il Y[X^^, = l,Y=l>Py 

i'=l j'=l 
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+ < n = = X n = i'^' = > (1 -Py) 

i'=l j'=l 

= Pn,,iPi„iPy + Pn„oPn„o(l " Py) ■ (68) 
Here, we denote the conditional default probabilities with the condition Y = as p^. g, 

Pn.,0 =< K.+i\ n X (1 - = 1 > • (69) 

n'=l ' 

Between p^. q and Pn-^o, the next relation holds. 

Pko = PkiPy + Pkoi'^ -Py)- (70) 
In addition, from the correlation between X^. and Y, we also have the next relations. 

Pn,,lPy = PifiPy + Pxy\jpifl{^ - Pifl)Py{^ - Py) (71) 

Putting these relations into ( l68l) . we can prove the next equations. 

i'=i j'=i 

= pUfiPifi + Piy X Pi2^\/K.,o(l -K„oX,o(l -Pn„o) X e-^("'+"^') (72) 

Using these relations, we calculate the inter-sector correlation and prove fl66l) . 

Next we need to check the validity of the above gluing process. We consider the 
two sector case K = 2 and their intra-sector parameters are set to be the same as A^, p 
and p in each sector. For the inter-sector default correlation pi2, we set pfy = ply = pu 
in the above multi-sector model. If the glueing process of the multi-sector model works 
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Figure 14. Comparison of Two sector MCB model: Directly Coupled Model vs 
Coupled by Y Model {py = 0.5). Plot of AP(n) = P2Ar(") - vs n/{2N)). 

N = 5, 10, 2Q,p = p = 0.03 and A = 0.3. Upper three curves have pxy = 0.01 and the 
bottom line with x symbol has pxy = 0.02. 



well, P2N{'n) should coincide with P/v+Ar(r2) of the coupled N + N MCB model in the 
previous subsection. Figure [TH shows AP{n) = P2Niji) — PN+N^n) with = 5, 10,20 
and A = 0.3. We set p = p = 0.03 and pi2 = 0.01. In addition we also plot the data for 
Pi2 = 0.02 and = 20. As the system size A^ becomes large, the discrepancy AP{n) 
increases. With the same system size A^ = 20, as the inter-sector correlation increases, 
the discrepancy also increases. As we have stated previously, these models obey the 
same conditions on the conditional default probabilities, however the default probability 
profile does not coincide. The glueing process by the auxiliary random variable Y may 
cause changes to the joint probabilities. We have not yet fully understood this point. 

With the present model, we study the effects of the inhomogeneous default 
correlation pmter 7^ Pintra on the default distribution function Pn{ti) and the loss rates 
D{i). We consider ideal portfolios which have the same default probability Pk = P 
and default correlation = p. The inter-sector default correlations are also set to be 
the same value as pij = pmter- Nk are also set to be the same Nk = Ng. We fix the 
total number of assets as A^ and compare P/v(t?.) and D{i) between the portfolios with 
different number of sectors K. Of course, between K and A'^, the relation N = Ns x K 
holds. We set Pinter < P- As the number K increases, the average default correlation is 
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governed by the inter-sector correlation Pinter and becomes weak. 

At the extreme hmit K = N case, each sector contains only one asset. In the 
inter-sector default correlation, Pn{ji) is given by the superposition of P/v(n|F = 1) and 
P/v(n|y = 0). The conditional default probabilities and ^^re given as 



Pi,o=P and Pl,=p + ^/p^■il-p). 



P (l-Py) (73) 



\1-P Py 



pIo=P- Vp^vP^^-^jY^^- (74) 

PN{n\Y = 1) and P/v(?7,|y = 0) are the binomial distributions 'Bm{N,pii) and 
Bin(iV,pto). 




Figure 15. Semi-log plot of PN{n) vs n. N = 100, K = 2,10,100, p = p = 0.03, 
Pinter = 0.01 and A = 0.3. Wc set Py — 0.5. Solid line represents K ^ 1 and p ~ 0.03 
case and the dotted line shows the curve for K ~ 1 and p = 0.01. 



Figure [I5] shows the semi- log plot of the default probability PAr(n) for K = 
1, 2, 10, 50, 100. We set the model parameters as p = p = 0.03, = 100, Pinter = 0.01 
and Py = 0.5. The solid curve plots the data for K = 1 and p = 0.03 and the dotted line 
shows the data for K = 1 and p = 0.01. As K increases, the data for each K departs 
from the solid line. At A' = 10, the data almost shrinks on the dotted line {K = 1 
and p = 0.01). The data for K = 50 almost coincide with those of A' = 100, whose 
PN{n) is given by the superposition of the two binomial distributions. This point is 
also a drawback of the present model. If the glueing process works perfectly, these data 
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should coincide with the homogeneous portfoho case K = 1 and p = 0.01. However, this 
discrepancy is inherent property of the modeL Contrary to the discrepancy in the two- 
sector case, the conditions on the conditional default probabilities p^^ i are different 
from those of the K = 1 homogeneous portfolio. They could cause the difference in 
PN{n). 




Figure 16. Plot of \AD{i)\ and D{i) vs i. AD{i) are between p = 0.01,0.05 and 
p = 0.03. N ^ 100, K = 2, 10, 100, p = 0.03, punter = 0.01 and A = 0.3. We set 
Py = 0.5. The solid curve shows the data of D{i) with p = 0.03. 



From the above discussions on PN{n) with different K, we think that the inter-sector 
default correlation Pinter is more important than the intra-sector default correlation p 
in cases of a large K. In the = 10 case, the PNin) are roughly given by those of 
the homogeneous portfolio with p = pmter- If one estimates the implied values of p and 
Pinter from the premium (or D{i)) of the portfolio with large K, this point is crucial. 
In Figure [ini we show \AD{i)\ for K = 10 between p = 0.01,0.05 and p = 0.03. We 
set the model parameters as in the previous figure. For comparison, we also plot D{i) 
for p = 0.03. D{i) represents the expected loss rate of the i-th tranche, the magnitude 
of |AD(i)| is important when one estimates the implied default correlation p from the 
premium of the tranche. We see that |AD(i)| with small i is small as compared with 
D{i). If we change p from p = 0.03, D{i) does not change significantly. It is difficult to 
obtain the implied values of p from the premium of the tranche with lower seniority. In 
contrast with medium values of i ~ 13, the magnitudes of D{i) and AD{i) are almost 
comparable. D{i) is sensitive to the change in p and it is not difficult to derive the 
implied value of p. 
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4. Implied Default Correlation 

In the last section, as a concrete example, we try to estimate tlie implied values of the 
default correlation from the premium of a synthetic CDO. We treat iTraxx-CJ (Series 2), 
which is an equally weighted portfolio of 50 CDSs on Japanese companies. The standard 
attachment points and detachment points are {0, 3%}, {3%, 6%}, {6%, 9%}, {9%, 12%} 
and {12%, 22%}. Table □ shows quotes on July 5,2005. The quote for the {0%, 3%} 
transhe shows the upfront payment (as a percent of principal) that must be paid in 
addition to 300 basis points per year. The other quotes for the other tranches are the 
annual payment rates in basis points per year. The index indicates the cost of entering 
into a CDS on all 50 companies underlying the index. The recovery rate R is 0.35. 

Table 1. Quotes for iTraxx-CJ Tranches on July 5, 2005. Quotes for the {0,3%} 
tranche are the percent of the principal that must be paid up front in addition to 300 
basis points per year. Quotes for other tranches and the index are in basis points. 
Source: Morgan Stanley Japan Securities Co. and Bloomberg 





{0,3%} 


{3%, 6%} 


{6%, 9%} 


{9%, 12%} 


{12%, 22%} 


Index 


5-year Quotes 


15.75 


113.25 


42.0 


30.5 


15.5 


24.55 



In order to get the implied default correlation for each tranche, it is necessary to 
relate the loss distribution function PN{n) to the premiums. The premiums are the 
present value of the expected cash flows. The calculation of this present value involves 
three terms |29]. We denote by Pfc(ai,aj^) the remaining notional for the {ai^au} 
tranche after k defaults. It is given as 

[ {an - aL)N k < KiV/(l - i?)] 

Pk{aL.aH) = \ aHN-k{l-R) K^^/(l - i?)] < A; < [a^^iV/(l - i?)] (75) 
i k> \aHN/{l - R))] 

Here, [x] means the smallest integer greater than x. For simplicity, we treat the 5-year 
as one period. The three terms are written as 

A = 5.0x < Pk{aL,aH) > e"^'"^ 

= 2.5 X {{an - aL)N- < Pk{aL, an) >)e-^t 

C = {{an - aL)N- < P^iaL, an) >)e-^i (76) 

where r is the risk-free rate of interest and we set r = 0.01. By considering the total 
value of the contract, one can see that the break even spread is given as C/{A + B). 

For the index, = 1.0 and = 0.0 and it is possible to estimate the average 
default probability p. Instead, we use the CDS data for each company and estimate the 
average default probability p and its dispersion AP^ as 

p = 1.8393%/5-year APa = 1.131%/5-year. 

With these parameters, the tranche correlations can be implied from the spreads 
quoted in the market for particular tranches. These correlations are known as tranche 
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correlations or compound correlations. As a pricing model, we use MCB, BBD and 
Gaussian copula models. For MCB, we use the following candidates. 

• Original MCB model (MCBl). N = 50,K =1 and A = 0.0. 

• Short tail MCB model (MCB2). = 50, K = 1 and A = 0.3. 

• Short tail MCB model (MCB3). = 50, AT = 1 and A = 0.6. 

• Disordered MCB model (MCB4). A^ = 25 + 25,A = 0.3, p^y = = py = p and 
Px = P + ^PdiPy = P — ^Pd- N + M MCB model with inhomogeneous default 
probability. 

• Two-sector MCB model (MCB5). K = 2,px = Py = p,Px = Py = P and p^y = 0.0. 
Assets are categorized in 2 sectors and pmter = 0.0. 



Table 2. Implied tranche eorrelation (%) for 5-year iTraxx-CJ on July 5,2005. 



Tranches 


MCB 1 


MCB 2 


MCB 3 


MCB 4 


MCB5 


BBD 


Gaussian 


{0%, 3%} 


11.79 


10.8 


9.96 


12.88 


21.2 


11.4 


13.8 


{3%, 6%} 


1.27 


1.18 


1.13 


1.36 


2.45 


1.26 


1.35 


{6%, 9%} 


3.16 


3.08 


3.09 


3.46 


6.32 


3.15 


3.23 


{9%, 12%} 


6.16 


5.95 


5.90 


6.65 


12.15 


6.11 


6.31 


{12%, 22%} 


9.78 


9.67 


9.90 


10.67 


19.97 


9.73 


9.46 



Table [2] shows implied tranche correlations for the 5-year quotes in Table [H We see 
a "Correlation Smile" , which is a typical behavior of implied correlations across portfolio 
tranches [T3] . 

We also find that the implied correlations are different among the models. For MCB, 
models with a larger A have a smaller correlation skew. The correlation for {0%, 3%} 
decreases with A and other correlations do not change significantly. If a probabilistic 
model describe the true default distribution, there should not exist any correlation skew. 
As a model approaches the true distribution, we can expect that the skew decreases. 
MCB3 is more faithful to the true default distribution. The skew of BBD is between 
MCB 1 and MCB 2, which is reasonable because the profile of BBD is between MCBl 
and MCB2 ( see figure [5]). Gaussian copula's skew range is larger than MCB 1, MCB 
2, MCB 3 and BBD. 

About the effect of APd, the implied correlations are considerably different between 
APd ^ (MCB 4) and APa = 0.0 (MCB 2). In the estimation of the implied correlation, 
we cannot neglect the fluctuation AP^. In particular, for the tranches {0%,3%}, the 
implied value is affected greatly by AP^. As has been discussed in [5], the increase in 
the dispersion of the default probability increases the loss in the equity tranche. It is 
necessary to increase the implied correlation to match with the market quote. 

In the case K = 2 (MCB5), the correlation level and correlation skew are very large. 
This means that if assets are categorized in many sectors and inter-sector correlation is 
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very weak, the loss distribution accumulates around the origin n = 0. In order to match 
with the market quote, a large intra-sector correlation is necessary. 

For the estimation of the model parameter A, we think that the resulting loss 
distribution should look closer to the implied loss distribution. That is, the range of 
tranche correlation skew should be small. Figure [T7] shows the tranche correlation vs 
A. As we increase A, the skew range becomes small. At A ~ 0.61, the range becomes 
minimal. We should calibrate A to be A = 0.61 ~ 0.62. 



Tranche Correlation vs k for 5-year iTraxx-CJ on July 5,2005 
0.12 r 1 1 1 1 1 1 1 r 
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Figure 17. Range pf implied tranche correlation for 5-year iTraxx-CJ on July 5, 2005 
vs A. 



5. Concluding Remarks 

In this paper, we generalize Moody's correlated binomial default distribution to the 
inhomogeneous portfolio cases. As the inhomogeneity, we consider the non-uniformity 
in the default probability p and in the default correlation p and Pinter- To treat the former 
case, we construct a coupled N + M MCB model and obtain the default probability 
function P/v-i-mI'^)- The inhomogeneity in p causes changes in the expected loss rates of 
the tranches with lower seniority. 

In order to treat the inhomogeneity in the default correlation, we construct a multi- 
sector MCB model by glueing multiple MCB models by an auxiliary random variable 
Y. We cannot take out the joining lines between the MCB models, for small portfolio 
and small Pinter, the construction works well. For the inhomogeneity in p, we divide 
a homogeneous portfolio into K sectors. We set pmter < P = Pintra and see the effect 
of the increase in K on Pj^{n). As the sector number K increases, the inter-sector 
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correlation Pmter becomes more important than the intra-sector default correlation p. 
With large K, the default correlation is governed by Pmter only. The CDOs, whose assets 
are categorized in many sectors, pmter should be treated more carefully than p = pmtra- 
In order to check the validity of the MCB model and our generalization method, 
more careful treatment and calibration should be done. We assume that p„ decays 
exponentially with n. With such a modification, the skew of the correlations diminishes, 
however, the skew remains significantly. Other models for p„ should be considered. For 
this purpose, it is necessary to study the implied loss distribution directly. Recently, 
Hull and White [281 129] developed a method to derive the implied loss distribution 
and to obtain the implied copula function from the market quotes of CDOs. The 
authors proposed a calibration method for p,„ from the implied loss function [23]. By 
incorporating this information in the MCB model's framework, we may have a "Perfect" 
correlated binomial default distribution model which reflects market quotes completely. 
We think that our generalization method provides important information based directly 
on the market quotes. 
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